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Abstract 



In this paper we describe gauge invariant multi-quark states generalising the 
A ■ path integral framework developed by Parrinello, Jona-Lasinio and Zwanziger to 

^ . amend the Faddeev- Popov approach. This allows us to produce states such that, in 

a limit which we call the ice-limit, fermions are dressed with glue exclusively from the 
fundamental modular region associated with Coulomb gauge. The limit can be taken 

K^ ' analytically without difficulties, avoiding the Gribov problem. This is illustrated by 

00 ' 

^j. ' an unambiguous construction of gauge invariant mesonic states for which we simulate 

^O ■ the static quark-antiquark potential. 

o 
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1 Introduction 

Yang-Mills theories are the cornerstone of the standard model. Their success is largely 
based upon perturbation theory where gauge fixing, implemented using the Faddeev- 
Popov method [1], plays a key role. However, Gribov [2] has pointed out that, at a 
non-perturbative level, the Faddeev-Popov method fails since there are always gauge equiv- 
alent (Gribov) copies which satisfy a chosen gauge condition. Initially this was shown for 
Coulomb gauge, but Singer has proven that this is in fact a general problem [3] . 

There have been various attempts to extend the Faddeev-Popov method to the non- 
perturbative regime. For example, it was suggested that the various Gribov copies, weighted 
by the Faddeev-Popov determinant, should contribute to the functional integral with al- 
ternating signs. This approach can be viewed as the insertion of a topological invariant 



into the partition function. Unfortunately, that topological invariant turns out to be zero 
in SU(A^c) Yang-Mills theory leaving us with the disastrous conclusion that the generalised 
Faddeev- Popov method results in physical observables being in indeterminate form [4-6]. 

This state of affairs is unfortunate as the need for non-perturbative gauge fixing is 
widely recognised as physically desirable. For example, Dyson-Schwinger equations are 
widely used in hadron phenomenology, and their construction relies on unambiguous gauge 
fixing, in particular in the infra-red regime. Using stochastic quantisation to by-pass the 
Gribov problem [7-9] , Zwanziger showed [10] that the tower of Dyson-Schwinger equations 
is unchanged but supplemented with additional constraints reflecting that gauge config- 
urations are confined to the first Gribov region. It turns out that the Green's functions 
solving the Dyson-Schwinger equations [11-13] appear to agree to a large extent with lat- 
tice simulations [14-16]. We note, however, that some initial discrepancies [17-19] in the 
infra-red [20] behaviour of Green's functions have been confirmed in large volume simula- 
tions [21,22]. It became clear only recently that these findings can be accommodated by the 
Gribov-Zwanziger approach when the Gribov-Zwanziger action is appropriately modified 
while preserving renormalisability and BRST invariance [23,24]. 

To go beyond the Faddeev- Popov method, we will here use an alternative construction of 
the partition function [25,26] which defines a gauge invariant action by integrating a weight 
function over the gauge orbit. This method has been studied on the lattice in the strong 
coupling expansion [27, 28] , and in numerical simulations of the weak-coupling regime [29] . 
The phase diagram was explored in [30] , and a phase transition from the weak to the strong 
gauge fixing regime was reported. Finally, it was argued in [31] that the gluon propagator 
displays gluon confinement. As we shall show, in a particular limit, which we call the 
ice-limit, the weight function constrains the gauge configurations to unique representatives 
of each gauge orbit. Altogether, these form what is called the fundamental modular region. 

In a Hamiltonian framework integration over the gauge group may be used to define 
projection operators onto the different non-Abelian charge (superselection) sectors of the 
Yang-Mills Hilbert space in the presence of external charges. This was first emphasised 
by Polyakov [32] and Susskind [33] and subsequently worked out in detail by a number of 
authors (see e.g. [34,35]). More recently Zarembo has used this approach to discuss the 
Yang-Mills mass gap, confinement and the interquark potential [36,37] (see also [38,39]). 
A thorough study of U(l) quantum mechanics along these lines may be found in [40]. 

In this paper, we introduce gauge invariant external fields (such as heavy quarks) 
through the projection techniques [32-37] into the above alternative construction of the 
partition function [25,26]. Using lattice regularisation, we will study the ice-limit where 
the fields are restricted to the fundamental modular region of Coulomb gauge. As an 
illustration, we will calculate the static heavy quark-antiquark potential. 



2 Non-perturbative gauge fixing 

2.1 The Gribov problem 

Recall that gauge fixing amounts to identifying the space 21/0 of gauge inequivalent config- 
urations (that is the space 21 of all configurations A modulo gauge transformations g E &) 
with a subset of the total configuration space: 21/0 C 21. The original idea of gauge fixing 
(due to Weyl, see [41] for the historical context) attempts at choosing a gauge 'slice' F 
(of configurations satisfying the gauge condition) to be identified with the physical config- 
uration space. While this works for the Abelian case, it fails for the non-Abelian theory 
due to the existence of residual gauge copies, as shown by Gribov [2], who was also the 
first to suggest a possible solution. As the copies only appear as one moves away from the 
perturbative small field regime and reaches what is called the 'Gribov horizon' it seems ap- 
propriate to just stay within its interior, i.e. within the Gribov region. Mathematically, this 
is defined as that neighbourhood of the classical vacuum {A = 0) where the Faddeev- Popov 
operator has a positive spectrum. It turns out, however, that this 'off-limits' prescription 
is not sufficient. Let us briefiy recapitulate the problem and its (formal) solution. 

Following 't Hooft [42] one may formulate the gauge fixing procedure in terms of dis- 
tance functionals S'fix[^^] = jl^^jl^ with || . ||^ an appropriate L^ norm. As shown by 
Semenov-Tyan-Shanskii and Franke [43] as well as dell'Antonio and Zwanziger [44] the 
extrema of S^x define the gauge condition while its Hessian (at the critical points) is the 
Faddeev-Popov operator such that the Gribov region is the domain of positive curvature 
containing A = 0, known to be convex and to cover all orbits [43,44]. Its boundary is 
the Gribov horizon where the lowest eigenvalue of the Faddeev-Popov operator vanishes. 
These authors also realised that there are copies remaining within the Gribov region, and 
one has to restrict configurations even further to the set A of global minima, 

A = {AeQi: Si,M] < SiiM% for all g E (3} . (1) 

Note that, by construction, this set is included in the Gribov region and hence the gauge 
slice. As pointed out by van Baal [45] one still requires suitable boundary identifications 
within dA endowing the subset A C 21 with the appropriate topology before it can finally 
be identified with the physical configuration space of gauge inequivalent configurations. 
In this context the latter is denoted the fundamental modular region (FMR), see [46,47] 
for reviews on this subject. We emphasise at this point that it is gauge invariant by 
construction. The details of the embedding 21/0 C 21, however, will depend on the gauge 
fixing (functional) chosen as its starting point. 



2.2 An alternative implementation of gauge fixing 

As we have discussed, a gauge fixing condition Xa[A] = can always be identified as a 
stationary point of a gauge fixing functional Sfixi^^] such that 

S^^[A3] = ^ Xa[A] = , (2) 



69''(x] 



where the fields 6°'{x) parameterise the gauge transformation, g{x) = exp{i^"(a;) t'*}, with 
t°- being the generator of the SU(A'"c) gauge group. The Faddeev-Popov approach is then 
based on the usual assumption [48, Chapter 16] that one can write 1 as 

Ve 6{xa[A^]) Det(^^^^ . (3) 

This, though, is not true non-perturbatively (see Appendix El for more details) as, in fact, 
the right hand side of ([3]) is zero. We will therefore use a different approach here [25,26], 
and, after reviewing it, we will study a series of examples. 

The starting point of this approach is the definition of a gauge invariant effective action 
5'cff[A] derived from the gauge fixing functional via the identity 

l = ^-SMA] fvge'^-^^'K (4) 

where T>g is the Haar measure on the gauge group. As it stands, this is a purely formal 
definition and one might ask if the right hand side of (jl|) is genuinely 1. In the continuum 
this question is hard to address, but using a lattice regulator it becomes clear that this 
really is a 1. To this end, we need to translate into a lattice formulation where the potential 
Afj^{x) is replaced by link variables U^{x) which transform under a gauge transformation 
(now conventionally written as VL{x)) according to 

U^{x) = n{x)U^{x)n\x + ae^) . (5) 

Inserting (jll) into the Yang-Mills partition function we obtain 

Z= JVU, Vn e^fl^t^"] e-^^«[^l e^--[^l . (6) 

For such a lattice regulated partition function, we may interchange the integration over 
the links [/^ and the gauge transformations fi: 



VnVU^ e^«''[^''l e-^^«[^l e^^^l^' = / VQVU^ e^«>=[^"l ^SMu"] ^s^mW^] 
ijvn\ /'l)t/^e^«''[^le"^^«[^le^™[^], (7) 



from the invariance of the action and the Haar measure. This means that we have been 
able to factor out the gauge redundancies into a volume factor in much the same way as 
the original Faddeev-Popov trick tried to do. However, as we shall see, this procedure is 
valid non-perturbatively. 

We will now investigate how this construction is used in three examples. 

2.3 Three examples 
2.3.1 The Christ-Lee Model 

The Christ-Lee partition function [49] is given by the two dimensional integral 

Z^,= j dFxe'c.i-\ (8) 

where S^^{x) is a function depending only on r = a/x^ + y'^. Gauge transformations are 
rotations through an angle 0, which we write x ^ x'^. Taking, for example, S^^{x) = —x^, 
one may check that Z^-^ = 27r x 1/2, where the 27r comes from the integral over the angle 
and 1/2 is the 'physical' partition function. We will consider the following gauge fixing 
functional 



S^^[x^] = --{x'^~vY , (9) 

where v is an external "gauge fixing" vector. The corresponding gauge condition 

dcj) (i ^ V X U , 



exhibits two Gribov copies: x"^ parallel or antiparallel to v. The FMR is given by those 
vectors x which (globally) maximise the gauge fixing action. In the present case, these are 
all vectors of arbitrary length parallel to v. Without loss of generality, we choose v = (1, 0) 
so that the gauge fixing condition becomes y = 0. The FMR is then given by the positive 
X-axis (with Gribov copies appearing on the negative x-axis). 
In analogy to (jlj), we find an effective action 

exp (SeffM) = fd(P exp {Si,,[x'^]) = 2nIoi^r) exp ( - '^ r' - '^Y (10) 

with Jo a modified Bessel function of the first kind. The effective action is manifestly gauge 
invariant as it depends only on r. We now insert the associated representation of unity. 



1 = exp ( - S,s[x]) fd(f) exp (^fi.[a^<^]) 
— - — - d(j) exp{Kx'^-v) , 



1 / ^^^' 
/ d(b exn ( k, x'^ ■ v) . 

27r. 



into the partition function ([H]) where it follows from the discussion in Section 12.21 that 

^CL = (/#) jd'x exp (^cl(^))^^^^^ , (12) 

using X ■ V = X. For S,^^ = — |icp, it is a straightforward matter to perform the d^x 
integral in flT^ and show that we recover the expected value of 1/2. The important point 
is that the Gribov problem does not hamper this calculation. This can be made explicit 
by considering the limit of large k, where the asymptotic behaviour of the modified Bessel 
function gives 



(fdA I dFxe'^-^^^J^e-'^^'-^K (13) 



Importantly, K,{r — x) > so that the k dependent terms of ( iT3l) are a Gaussian regulari- 
sation of the delta function. The support of the delta function arising in the k ^ oo limit 
are those vectors x for which 

n{r — x) = =^ y = and x > 0. 

The condition 2/ = corresponds to our chosen gauge condition, but the condition x > 
restricts us to only the FMR, i.e. the Gribov copy at x < is not seen by the partition 
function. Explicitly: 

lim J— e-''^"-^) = |x| ^(x) 5{y) , (14) 

K— »oo y 27r 



where is the Heaviside step function. Using this in flT3|) our final result for the k ^ oo 
limit is: 

Z^^ = ( fd(yj fd^x e^cLW^(3.) |a;| s{y)^ (15) 

Our approach has not only correctly produced the gauge fixing constraint in terms of 
the 6 function and the Faddeev-Popov determinant |x|, but also the correct "horizon 
function" [50,51] ^(x), which singles out the FMR to the right of the Gribov horizon at 
X = 0. In Figure[T]we give contour plots of the k dependent fraction in flT^ . at k = 16. We 
clearly see the FMR emerging as the domain of support. The discussion of a more general 
class of gauge fixing functions is left to Appendix [Bl 

2.3.2 Example: U(l) Landau gauge 

As a second example we consider U(l) gauge theory in Landau gauge. Although this does 
not have a traditional Gribov problem, zero-modes must be eliminated in order for the 
Faddeev-Popov determinant to be non-zero. 

The gauge fields A^{x) change under gauge transformations as 

fi(x) = e^^(^) , A)^(x) = A^{x) + d,\{x) . (16) 
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Figure 1: Contour plots of the k dependent fraction in ( IT2l) . at k = 16. The support of 
this combination of Ses and S^x, and thus of the partition function Z^^, is restricted to the 
FMR at large k. 



We choose Dirichlet boundary conditions: 

Af,{x)\xedv = 0, X{x)\x(^ov = Ao 

with Ao being constant. The gauge fixing action for Landau gauge is given by 



^fix[A] = -m^ fd^x A^A^ix) , 



(17) 



:i8) 



where the "mass" m acts as a gauge fixing parameter. This gauge fixing action generates 
the Landau gauge condition via 
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S\{x 



-Si^^A' 







d^Af, = 



Given the boundary conditions flTTl) . it follows that any A(a;) may be written 

\(x) = Ao + A(a;) . 



(19) 



(20) 



where Aq is the surface value of A(a;) and A(a;) vanishes on dV (A may be decomposed as a 
sum over the non-zero Fourier modes of the Laplace operator). The measure on the algebra 
elements A is inherited from the group, so that VQ = dAoI^A. Note that the gauge fixing 
action ( TT8l) is invariant under constant gauge transformations since A ^ A + dX = A + dX. 
We must adopt some prescription to deal with the zero modes and to this end we will 
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include in our measure a delta function S{Xq) which kills the zero mode and makes the 
d'Alembertian D invertible. We are therefore led to propose the following form for the 
effective action: 

^SMA] ^ f^j gxp (5fix[#]) = det-^/\-m^n) exp ( - Si,^[A^]) . (21) 

The transverse gauge field A^ is defined by 

Al=(^V,,-^d,d}jA'^, (22) 

and is gauge invariant. This result may be used in ([7]) to obtain the partition function of 
U(l) gauge theory. We evaluate 

^fix[A] - S,s[A] = -m" fd'x A'^A'^^ + ^ logdet(-m2n) , (23) 

where the longitudinal part of the field is defined by A^ = A — A^ . Inserting our repre- 
sentation of unity into the partition function we find that for this U(l) theory 

Z = det-^/^(-m2n) fvA^ exp ( - Sym[A] - m^ fd^x A^A^^) . (24) 

Note that the gauge fixing parameter m acts as a mass for the longitudinal gauge fields. In 
the large mass limit, these decouple from the partition function leaving us with transverse 
fields only. 

2.3.3 Example: SU(2) and weak gauge fixing 

In our final example we consider SU(2) Yang-Mills theory, in lattice regularisation, with 
the gauge fixing functional 

5fi.[f/^] = ft:^trf/^^(x), fi=l..A, (25) 



X,fl 



which implies (lattice) Landau gauge upon extremisation. When applied to the partition 
function Z, the (local) maxima of Sfi^ give the dominant contributions at large k. We refer 
to this as the 'strong gauge fixing' limit. 

As we have illustrated with our previous models the use of the effective action approach 
is not limited to the large k, regime. Here, we take k <^ 1, the 'weak gauge fixing' limit, 
and calculate the effective action perturbatively in k. Expanding the defining equation (jl]) 
and noting that S'gx is of order k, we find, up to fourth order in k, 

So^U] = l{Sl) + ^[{SL)-S{Siy]+0{K'), (26) 

(S^^) := [VQS^U% (27) 



For SU(2), terms with an odd power of k vanish upon integration over Q. The term 
quadratic in k is independent of the hnks since there is no gauge invariant combination of 
two hnk variables apart from a constant: 

(SL) = ^N, (28) 

where A'^; is the number of hnks on the lattice. The calculation of the fourth order term is 
tedious but straightforward. We find: 

(SL) - 3 {SLY = y^^ + ^ E l''Pr>[U], (29) 

p 



where the sum extends over all plaquettes, Pp. Inserting (|29|) and (l2g|) in fizol) . we finally 
obtain: 

"1 /K\2 1 /K^ 



Ses[U]=Ni 



2 



©^-l^(f)>Kf)*Ei-.l-H-(A (30, 

For sufficiently small k, the effect of S'cff[f/] is just to correct the coefficient jS of the 
plaquette in the Wilson action, showing explicitly that ^efiff/] is gauge invariant. Terms 
such as the 1x2 Wilson loop appear at order k^. 



3 Non-perturbative dressing 

In this section we will introduce gauge invariant heavy quarks into the above approach 
to the partition function. We will adopt a Schrodinger representation, considering gauge 
invariant states constructed in a single time slice. Transition amplitudes between such 
states will be discussed in section HJ We begin by briefly reviewing the construction of 
gauge invariant charges [52]. 

A gauge invariant charged state | Q ) may be constructed from a fermionic state q{x) \ ) 
by 'dressing' the latter with an appropriate function h[A] of the gauge field A, 

\Q):=h[A]{x)q{x)\q). (31) 

From the transformation properties of the fermion, q{x) -^ Q{x)q{x), requiring gauge 
invariance of the state | Q ) implies that the dressing h[A] transforms as 

h^[A]{x) = h[A]{x) Q\x) . (32) 

Dressings may be constructed through a field-dependent gauge transformation which ro- 
tates a given field A into a gauge x[^] = . The dressing factor is not unique and should 
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be selected according to its physical properties [53,54]. For example, the dressing describ- 
ing a single static quark in SU(A^c) comes from the field dependent rotation into Coulomb 
gauge defined by 

di(hAih-^ + ihdih-A =0 . (33) 



From this example it is clear that the above dressing approach is sensitive to the Gribov 
problem of Coulomb gauge: although a gauge invariant charge may be constructed pertur- 
batively from a solution of ( l33l) . Gribov copies imply that the dressing factor is not well 
defined non-perturbatively [55]. This offers a route to understanding confinement as this 
lack of single observable quarks is in agreement with experiment. The open dynamical 
question is how the Gribov ambiguity produces the physical scale of hadronic multi-quark 
systems. 

3.1 Dressing by projection 

In this section we review the projection, or group integration, method of [32-37] to con- 
struct gauge invariant states, and then make the connection with the partition function 
above. Consider, for illustration, the group integral 

h[A]{x) := fvn n{x) e^[^"] , (34) 

with an, a priori arbitrary, weight function exp(Vr[A]). It may be checked that (!M|) obeys 
the transformation law ( |32l) and therefore gives a dressing for a single charge. Dressings for 
multi-fermion states may be similarly constructed by inserting the appropriate factors of 
Q or fi"'' under the group integral. In this paper we propose the weight functional W = ^gx, 
which will allow us to combine the projection approach with the construction of the gauge 
fixed partition function discussed in section 2. 

We have made two assumptions in writing down these dressings: (i) the group inte- 
gration exists (this is certainly the case in, e.g., lattice gauge theory), and (ii) the group 
integration does not yield a vanishing result for h[A]. 

The latter assumption is more restrictive. Indeed, we will see below that (ii) is not 
fulfilled for the Coulomb dressing, and that (1341) needs modifications for this important 
case. In fact we will see that the required changes allow us, importantly, to attribute a 
global charge to our locally gauge invariant states. 

Let us first look at an Abelian example of this. 

3.2 Example: U(l) dressing and zero modes 

In this example we construct a U(l) dressing from the Coulomb gauge fixing functional. 
To illustrate the importance of flat directions, i.e. invariance under some class of gauge 
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transformations, in S'gx, and how they are related to global charge, we will work in a finite 
cubic volume V = L^, with boundary dV (recall our states are defined in an initial time 
slice). We again take Dirichlet boundary conditions on Ai{x), and A(a;) constant on dV 
so that we may write A(a:;) = Aq + X{x), analogously to (l20ll . 
The gauge fixing functional for Coulomb gauge is 

Sa^lA] = -m fd'x Aiix)Mx) . (35) 

V 

We would now like to perform the U(l) analogue of the group integration ( IMl) to construct 
a dressing for a U(l) fermion. Again, the measure is VQ = dX^VX. However, S'fix[^] is 
invariant under constant transformations. It therefore does not see the zero mode Aq, as 
'S'fixi^'^] = 5'fix[^'^]- It follows that the integral over Aq vanishes, and therefore so does the 
dressing. To avoid this we once again add 6{Xo) to the group measure and find 

h[A]{x) = f dXoVX 6{Xo) exp {iX{x) + Sfi^[A^]) 

= fvX exp ( tX{x) + Sfix[A^]) ^^^^ 

= AT exp (5fix[A^]) exp {tA~'diA{x)) . 

Here J\f = det^ (— -mV^) exp [A^J^/4m] is a normalisation factor which is IR finite but 
UV divergent. The inverse Laplacian A"^ is well defined on Ai as the fields obey Dirichlet 
boundary conditions. The transverse field Af is gauge invariant while the longitudinal 
term A~^diAi is non-invariant. One may check using any of the expressions in ( l36l) that 
the behaviour of h under a gauge transformation a{x) = oq + ci{x) is 

h[A'']{x) = h[A\{x) e-'^'^'^K (37) 



We see that excluding the zero modes from the integral flHbj) has given us a dressing which 
transforms as in ( l32l) under local transformations, but is insensitive to global transforma- 
tions. It follows that the dressed fermion state transforms with a constant phase, 

\Q») = h[A] (a;)e-*"("') e'"^''k\ ) = e^"° \Q) , (38) 

which allows us to assign a global charge of one to these states. We now consider the 
analogous Coulomb dressed state in SU(A'c)- 

3.3 Example: multi-quark states in SU(A^c) 

The Coulomb gauge fixing functional in SU(A^c)5 using lattice regularisation, is 

5fix[f/^] = «: Yl ]^trf/,^(x). (39) 

x,l=1...3 ^ 

12 



As in the U(l) example above, S'gx is invariant under constant gauge transformations C 
because of the cychc invariance of the colour trace. The corresponding dressing obeys 

h[U] = I ViSlC) n e^fl>=[^^"''^l = I ViSlC) iSlC) e^fl>=[^"''l C^ = h[U] C^ , (40) 

and therefore vanishes. The difference between this and the U(l) example is that the 
measure on SU(A^c) does not allow us to separate out the constant gauge transformations 
and preserve the dressing property (1321) . For the remainder of the paper we confine ourselves 
to multi-quark singlet states which are invariant under global transformations. Let us 
concentrate on physically relevant states such as mesons and baryons adapting the methods 
of [35] for our purposes. 

For example, the dressing for a baryonic state in SU(3) would be 

h^^'^UUM^y.x) = y ^^ e,,tfi„(z) n^M ntn.{x) e^«''[^"l . (41) 

The baryonic trial state is invariant under global transformations C since 

The dressing for a quark-antiquark state appears as a straightforward generalisation of 

dMD, 

.-I^|(,..)^/Z)«..(.)n(.)e-'-'. (42) 

This dressing transforms homogeneously, 

h('\U'']{y,x) = G{y)h^'\U]iy,x)G^{x) , (43) 

and is invariant under constant transformations. The quark-antiquark state 

\QQ):=q{y)h^'^[U]{y,x)q{x)\0) (44) 

is therefore gauge invariant. 

In the strong gauge fixing limit, the dominant contribution to the dressing comes from 
the gauge transformation which globally maximises the gauge fixing functional. This is the 
transformation fiFMR[^](a^) which transforms a given A to its gauge equivalent configura- 
tion in the FMR, so for the mesonic states 



\QQ) ~ ?(?/) ^fmr(?/) ^FMR(a;) q{x) |0) e 



SrAA] 



(45) 

K>1 



In this case, the external quarks are dressed by those gauge transformations which rotate 
the gauge field into the FMR. In the following section we use this result to derive our main 
finding - we will be able to carry out the strong gauge fixing limit analytically without 
constructing the FMR explicitly. 
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4 Applications 

4.1 The heavy quark potential 

We now use the formalism above to address, in SU(2) gauge theory, the static potential 
of a heavy quark-antiquark pair made gauge invariant via Coulomb dressing. The ground 
state energy in this channel depends on the quark-antiquark separation r and equals the 
static potential V{r). Using the "mesonic dressing" (H2l) . we calculate the matrix element 
(where r = \x — y\) 

P^^^^T) = j-^^^^^^ , (46) 

\QQ). = qiy)h^'\y;x)qix)\0), 

where H is the Yang-Mills Hamiltonian. Note the dependence of the states on the gauge 
fixing parameter k. For Coulomb dressing, the gauge fixing functional is separately defined 
at each time slice and depends on time t through the time dependence of the background 
field U^ix): 



5fi.[f/](t) = «:^itrf/Ka;,t). (47) 

X, I 

The static quark-antiquark potential V{r) can be extracted from the large T limit of p(r, T) 



via 



p.(r,T)^|(2|QQ)J%-^('-)^, (4J 



where |2) is the true ground state in the quark-antiquark channel. Our aim will be to 
consider the strong gauge fixing limit /t — > oo for which the quark-antiquark trial state is 
dressed with glue from the FMR (see (l45l) ). Given that the static heavy quark propagator 
from time to T is proportional to the so-called "short" Polyakov line [56], i.e., 

P[U]ix,0,T)= II Uo{x,t). (49) 

te[o,T] 

we finally obtain for p^ in (H6l) 

p,{r,T) = (n\x,0)P[U]ix,0,T)n{x,T)n\y,T)P^[U]iy,0,T)n{y,0)')^, (50) 

where the latter expectation value is defined by 

(0)p = N^' I VU^ Vn O eSYMlu]+s^iu»] ^ ^ ^ f j^jj^ ^SymIu]^ (5^) 

Here the integral VQ is initially taken over the gauge transformations in the initial and 
final time slices, defining the states. However, with a properly normalised Haar measure 
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this may be extended to the entire lattice, which is how we are to understand VQ here and 
below. We have also introduced the total gauge fixing action by 

^F[t/] = 5fix[t/](0) + Sfi.[f/](T). (52) 

4.2 The weak gauge fixing limit 



The gauge fixing action flozl) contains the gauge fixing parameter k which controls the 
overlap of the Coulomb dressed quark-ant iquark state with the true ground state. In fact, 
for K = 0, our ansatz for the trial state vanishes, leaving us with a null result for pi^{r,t) 
in ( H6|) . For small but non- vanishing k, the overlap is non-zero, and we should be able to 
recover the full static potential in this limit of weak gauge fixing. 

Let us consider pK.{r, T) in fl50p in leading order of an expansion with respect to n. For 
this purpose, we will first perform the integration of Vt using techniques familiar from the 
strong coupling expansion of gauge theories. To leading order, we may expand the gauge 
fixing action: 

e^^[^"l = n(l + it^f^O' xG(a;,0), (a.,T), / = 1...3. (53) 

X, I 

Subsequently, we perform the integration over the gauge transformations restricting to the 
minimal power of k for which 053p produces a non- vanishing result. This technique is 
standard textbook material [57] so that we only quote our final result, 

/ K\ 2(r/a+l) 

p,{r,T) = [-) W{r,t), (54) 

where W{r,t) is the rectangular Wilson loop of spatial and temporal extent r and T, 
respectively. The static potential can be recovered in the standard fashion 

W{r,T) ^ |(2|a)|2e-^(^)^, (55) 

where \a) is the axially dressed quark-antiquark trial state [58,59], usually associated with 
a chromo-electric string joining the quark and antiquark. 

There are two important implications of (15^ : (i) the k dependent prefactor does not 
infiuence the static potential V{r) implying that we analytically find the correct potential 
for small but non-zero k; (ii) the overlap of our gauge invariant trial state with the true 
ground state is quite poor, at least in the weak gauge fixing limit: 

/K\2(r/a+l) 

K2|gg).r = y mo)\\ m 

The overlap is even lower than that associated with "thin" Wilson lines. For practical 
calculations using lattice gauge theory, we will seek values k, of order one. 
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4.3 Strong gauge fixing and the ice-limit 

In the last subsection we have shown that we can recover the static quark-antiquark 
potential in the weak gauge fixing limit without being hindered by the Gribov problem. 
We now show that this also holds for strong gauge fixing. To contrast our approach with 
the common difficulties previously encountered we will first show how the Gribov ambiguity 
reappears for large values of k. However, we will then demonstrate how to circumvent the 
Gribov problem in the dressing approach by taking the limit k -^ oo in an appropriate 
and novel fashion. 

Consider the matrix element /^^(r, T) given by the functional integral (l50l) . and how we 
might attempt to perform the orbit integrations VQ before the integration over link vari- 
ables, Ufj^. In other words, we integrate over Q treating the links t/^ as a fixed background 
configuration. Hence, in this step, the matrices fl are considered as the fundamental de- 
grees of freedom of a theory with probability weight exp{S'F[t/^]}. For SU(2) gauge theory 
we have 

Q = ujq + iuj ■ T , ul + UJ ■ UJ = 1 , 

identifying these degrees of freedom as 4-dimensional vectors (spins) of unit length such 
that the associated partition function possesses a global 0(4) symmetry. These spins 
interact only with their nearest neighbours. As the background links provide a nontrivial 
"metric" for these interactions the partition function describes a spin-glass. In the strong 
gauge fixing limit, k, —>■ oo, this spin-glass will approach its ground state. However, for 
generic background fields, there exists a variety of highly degenerate near ground-states 
leading to frustration of the system. Finding the global maximum of S'f[?7^], 

Sf[U^] ^ max. (57) 

amounts to identifying the true ground state of this spin-glass. It is well known that 
dealing with this problem is extremely costly and beyond the scope of standard numerical 
techniques such as importance sampling. Hence, the Gribov problem has reappeared in 
disguise, as the problem of simulating a spin glass at low temperatures. 

The crucial idea to avoid any Gribov (or spin-glass) problem is to trivialise the orbit 
integration over Vt by factoring it out altogether. This can be done as follows. First note 
that in the matrix element p^ from fl3Ul) we may write 

fi"f(a;,0)P[f/](x,0,T)fi(x,T) = P[[/^](a;, 0,T). 
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t-T 




Figure 2: Illustration of the ice-limit: action density and orientation of the finite length 
Polyakov lines (vertical lines). The low action planes (at t = and t = T) are clearly 
visible. 

Hence, using the gauge invariance of action and Haar measure, we obtain for iVp^ 

VU^, Vn n^{x, 0)P[U]{x, O, T)n{x, T) n^{y, T)P^[U]{y, 0, T)n{y, 0) e^YM[c/]+5F[c/"] 
VUf, Vn P[U^]{x,0,T)P^[U^]{y,0,T) e^™[^]+^F[^"l 
VU^ Vn P[U^]{x,0,T)P^[U^]{y,0,T) e^YM[t/"]+Sp[[/"] 
= (f'Dn) /"pf/^P[[/](a;,0,r)pt[;7](^,0,T)e^™[^l+^p[^], 

and, indeed, the gauge group integration factors out as a trivial volume factor, J Vfl = 1. 
Thus, our final answer for p^ becomes 



p,{r,T) = N~' jvU,P[U]{x,0,T)P^[U]{y,0,T)e 



Sym[U]+St,[U] 



N 



Pt/^e^™t^l. 



(58) 
(59) 
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Figure 3: The static quark-antiquark potential from simulations with k = 2.0 (left) and 
in the ice-limit k = oo (right). 

It remains to discuss the strong gauge fixing limit. In contrast to ( 1571) we now need to find 
the maximum of the gauge fixing action with respect to the links, f/^, 



Sy[U] 



max. 



(60) 



Crucially, this does not constitute a spin-glass problem as the global maximum of ( !60|) is 
easily obtained: 



Ui{x,t = Q) = l, 



Ui{x,t = T) = l, 



1 = 1... 3. 



(61) 



Hence, for infinitely large k, the spatial links of the time-slices t = and t = T are both 
frozen to the perturbative vacuum. We therefore call the limit /t — * oo the "ice-limit" of 
the integral dressing approach. The (finite length) Polyakov lines start and end on the 
frozen time slices. Figure [2] shows the action density in the cube consisting of the time axis 
and two spatial directions. 

We now present our numerical results. For finite values of k, (with n not too large 
to avoid ergodicity problems in the spin-glass limit), we generated configurations {U^,Q) 
corresponding to the partition function 



X)U qSym[U]+Sp[U] 



r'f/„r'fie^™[^l+^^[^' 



Sym[U] = PY. l^^P^Mi^) 



(62) 
(63) 



x,^>u 
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recalling f VQ = 1. In order to ensure that the excited states are sufficiently suppressed, 
the Euclidean time T must be chosen large enough. Here, we have considered all values 
T > Tth for the straight-line fit of 

-lnp{r,t)^ = V{r)T - ln|(2|QQ)|2 (64) 

where the asymptotic form (1481) of Pk(^) t) was used. Our final result for V{r)a using a 16^ 
lattice, P = 2.3 and 10000 independent configurations is shown in Figure [H] (left panel). 
The line is a ffi of the Tth = 4 data to 

V(r)a = aa^--a/(r/a) + Vo, aa^ip = 2.3) = O.U , (65) 

a 

where the known value for the string tension in lattice units was used. We find that Tth = 4 
is sufficient to decouple the excited states. A good agreement with the potential obtained 
with standard overlap enhancing techniques is observed. 

We finally study the ice-limit k — ;> oo. In this case, configurations were generated using 
the standard partition function except that the spatial links in time slices t = and t = T 
were fixed to unity. Using f3 = 2.3 and a 16^ lattice, the static potential was extracted from 
10,000 independent configurations in the ice-limit. The result is shown in Figure [S] (right 
panel). We observe that rather large values for Tth are required to decouple the excited 
states in this case. One needs Tth = 6 to achieve good results for V{r). Note, however, that 
no smearing was involved and gauge fixing ambiguities ('Gribov noise' [56]) are absent. 

5 Conclusions 

Non-perturbative gauge fixing is a key ingredient of many approaches to Yang-Mills theory. 
In order to have full analytic or numerical control all such approaches must confront and 
understand the Gribov problem as this unavoidably arises in any direct imposition of a 
gauge fixing condition. In this paper we have seen how to by-pass this problem both in 
the definition of an unambiguous partition function and in the construction of suitable 
mesonic states. To this end we have designed a novel and unified framework that combines 
the properly gauge fixed path integral introduced in [25, 26] and a generalised dressing 
approach based on group integration to construct gauge invariant states. The feasibility 
of this framework has been explicitly demonstrated for various examples ranging from the 
Christ-Lee model, U(l) gauge theory to SU(2) Yang-Mills theory. 

From our new vantage point, the emergence of the Gribov problem could be traced 
back technically to a problem associated with the order of two integrations extending over 
the gauge orbits and field configurations, respectively. If the gauge group integration is 
performed before the average over the gluon fields, the Gribov problem arises as the problem 
to find the ground state of a spin-glass. However, since gauge invariance is manifest in 
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our approach, the Gribov problem can be avoided by interchanging the integrations over 
gauge group and gluon fields. In this case, the ('ice') limit of strong gauge fixing can be 
performed analytically, and the external matter fields become properly dressed by gauge 
transformations unambiguously connected to the FMR of Coulomb gauge. The numerical 
feasibility of the method was finally demonstrated by a lattice calculation of the static 
quark antiquark potential from trial states in the ice limit, living on (initial and final) time 
slices frozen to the perturbative vacuum. 

Acknowledgments: The authors thank Andreas Wipf for helpful discussions. The 
numerical calculations in this paper were carried out on the HPC and PlymGrid facilities 
at the University of Plymouth. 

A Faddeev-Popov method and Gribov problem 

The standard approach to (nonabelian) gauge fixing is the Faddeev-Popov method. Let 
us briefiy recall it in the slightly generalised framework of [4-6]. In this modified approach 
one generalises the usual representation of unity ([3]) to the topological invariant 



N[A -- 


J ^ ^ \ 5e''{x) 


M[A]abix,y) -- 


5'5fix[AT 

5e<'{x)5e\y)' 



(66) 

(67) 

Under the idealising assumption that there is a unique solution to the gauge fixing condition 
(featuring in the 5-function in (!66l) ) the standard identity ([3]) is recovered, 

N[A] = 1 , (68) 

with a gauge invariant Faddeev-Popov determinant, DetM[y4] = DetM[y4^]. This may in 
turn be used to remove the gauge group volume from the partition function, 

Z = JVA, e^[^] = JVA, Ve' DetM[A] 5(^^^)e^--[^] 

exploiting the invariance of the Haar measure, the action and the Faddeev-Popov determi- 
nant. Interchanging the order of integration and renaming A^ -^ A^ finally yields 



[5S^AA] 



Z=(^jve'^ JvA,BetM[A]6[^: 



(x) 
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qSym[A]_ (gg^) 



The latter equation is the desired result: the trivial factor from the gauge degeneracy 
has been factored out from the partition function, and the residual integration can be 
straightforwardly evaluated by means of perturbation theory. 

The crucial observation, due to Gribov [2], is that the gauge fixing condition has many 
solutions. Hence, the group integration in (l66ll becomes a sum over all residual gauge 
transformations which cast a given background field ^^ into one of its Gribov copies. Due 
to the compactness of the group integration this implies that flUHj) is actually replaced by 



N[A] = , (70) 

so that also this generalised Faddeev-Popov approach remains ill-defined [5, 6] . 

B Christ-Lee model revisited 

In order to make contact with the Faddeev-Popov method we need to evaluate, from (fT2l) . 
the expression 



Scs[x] = Sfi^[x] -In / , 



gSfl,[cc"* 



for large values of k, and for an arbitrary Sa^. To this end, let (po denote the gauge 
transformation which transports the vector x along its orbit to the global maximum of the 
gauge fixing action. For large k, we may use a semi-classical approximation to evaluate Ses, 
via: 



In f d<P e^fl"!^^] = 5fix[x'^«] - - In M[x\ + In v^, (71) 



where M is the gauge invariant "Faddeev-Popov matrix". It follows that ( 1121) becomes 

Z = (27r) [ (fx e^CL(^) gSfl,[cc]-Sfl,[3;^0] ^^lnM[x]-lnV2^ _ (^-^3^ 

For S'fixia^'^] as in (^ with v = (1,0) and for a given vector x = (x, y), (t>o[x] is defined by 



X cos 4)0 — y sin 0o -^ max. 



This implies 



X cos(j)o[x] - y sin 00 W = r, Sii^[x'^°] = --{r - if. 



and indeed the Faddeev-Popov matrix is gauge invariant, M[x] = nr. 
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Returning to a general 3^^, consider the weight factor P{x) = exp{Sfi^[x] — S{i:^[x'l'°]} . 
This will restrict the field (here, x) integration to the FMR for large values of n. To show 
this, it is convenient to decompose the vector x into a part x^^^^^ G FMR and a fiuctuation 
along the gauge orbit. 



x = x'P 



FMR 

In contrast to the Faddeev-Popov approach, this must not be done for the 2-dimensional 
space as a whole but only locally for x close to the FMR. Since 



dSfiJx^ 



0, 



dip 
we are led to 

V dip ^ecc^Mi^J 
The weight function h{Xj,^^^) can be calculated by integrating ip over a small interval 
around zero. We will assume that the map x ^->- {x^^^ , ip) is invertible for cc in a region 
around the FMR. We stress, however, that this depends on the gauge choice and might not 
be the case for more general settings thus hinting at a shortcoming of the Faddeev-Popov 
approach. Under the above assumption, we use 

d^M{x^^^)5\ — ■^—- 



tcGtCp 



and find 

h{x,^^) = [dip Mix,,,J e^-[-l-^-[-^''] = Mix,,J I dip e-^*^(-™a)^^ = ^2nMix, 



(75) 



Inserting ( |75l) and (ffll) in (ff3l) . we finally obtain the desired result, 

Z = (2.) / d'x e^c.(^) M(a._) 6 (^%M ) , (76) 

which is the starting point of the Faddeev-Popov method. 
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